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, Abstract. We study unfrustrated spin Hamiltonians that consist of commuting tensor 

T-H ' products of Pauli matrices. Assuming translation-invariance, we observe that the Hamihon- 

, ian is described by a map between modules over the translation-group algebra, so homologi- 

' cal methods are applicable. We show that in any dimension every point-like charge appears 

as a vertex of a fractal operator, and can be isolated with energy barrier at most logarithmic 
^ ' in the separation distance. For a topologically ordered system in three dimensions, we prove 

that there must exist a point-like nontrivial charge. We also develop tools to compute the 
' ground state degeneracy and to handle local unitary transformations. 
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Commuting Pauli Hamiltonians form a small class of Hamiltonians that are consisted of 
products of Pauli matrices such that each term commutes with any other terms. Classical 
examples are the Ising models in one or two dimensions. Albeit its simplicity of the energy 
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spectrum, there are many intriguing models in this class for which the long range entangle- 
ment of the ground state plays a very important role. Prototypical is the Kitaev's toric code 
model [1], which has been a solid testbed of ideas for topologically ordered systems. 

The topological ordered models exhibit, as the name suggests, many properties that are 
insensitive local changes or defects. They had been discussed for the states of the frac- 
tional quantum Hall effects and the spin liquids; see e.g. Wen [2]. Perhaps, the most 
well-defined characteristic of the topological order is the local indistinguishability of the de- 
generate ground states; two different ground states gives the same expectation value for any 
local observables. (Note that this characteristic is not directly applicable to e.g. the topo- 
logical insulators [3], for which certain symmetry properties distinguish them from trivial 
phases.) Due to the local indistinguishability, the topologically ordered systems are thought 
to be candidate media on which quantum information processing is performed. As a special 
application, the topologically ordered system can be used as a quantum memory, just like 
the ferromagnetic system is used as a classical memory. 

However, the quantum memories in the topologically ordered systems often suffer from 
thermal instability. For example, the toric code model has point-like excitations, which can 
freely propagate by external noise from the thermal bath. Although a local operator can 
never access to the ground space, their accumulation may be able to. Indeed, by the thermal 
fluctuation, a ground state is often mapped to a different state, and the anticipated protection 
of the stored quantum information is not viable. The excitations that affects the stability of 
the quantum memory may be called "topological charges". A charge is an excitation that 
cannot be created alone locally but can be created with some other excitations. Indeed, the 
4D toric code |1] has no charge at all, and can be used as a quantum memory whose failure 
probabihty decreases exponentially with the system size at low enough temperatures [5l [6] . 

The situation in three dimensions is more subtle but interesting. Models like 3D toric code 
model have charges that can freely propagate across the system by the interaction with the 
thermal bath, thereby two different ground states become mixed. On the other hand, as in 
the cubic code [Zl [HI E] , there can exist charges that cannot propagate by any means. This 
class of models provides modest reliability as a quantum memory at nonzero temperature. 
However, the scaling of the memory time, until which the system is reliable as a memory 
medium, is not as favorable as it is for the 4D toric code model; the memory time grows 
with the system size according to a power law whose exponent is proportional to the inverse 
temperature, provided that the system size does not exceed some critical value determined 
by the temperature. 

The very existence of the charges seems to adversely affect the memory time. In order to 
have a quantum memory, one needs to devise a read-out procedure explicitly — a classical 
analog is the measurement of the average magnetization of 2D Ising model. Though the 
charges may not propagate, they can be separated arbitrary far from their partner charges 
at a modest energy cost [8j. It may sound contradictory, but crucial is that a set of charges 
can expand only in a highly restrictive way. No good read-out procedure (sometimes called 
decoder) is known for the configurations of far separated charges. This should be contrasted 
to the 4D toric code model, in which any excited state consists of several loops. The large 
loops are suppressed by the Boltzmann factor, and the small loops can be almost perfectly 
treated by the read-out procedure. In short, the large entropy due to the point-likeness of 
the charges would likely drag the ground state to a hard-to-decipher state. See the discussion 
in 0. 
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Apart from the issue of the thermal stabihty and the possibihty of quantum memory, 
the cubic code model apparently necessitates new tools to analyze it. When defined on a 
finite system with periodic boundary conditions, it shows exotic dependence of the ground 
state degeneracy on the system size. The degeneracy is sensitive to the number theoretic 
property of the linear system size L. For example, when L is a power of 2, the degeneracy 
grows exponentially with L, but becomes a constant if L = 2^ + 1. This was a numerical 
observation, and was not rigorously treated 

Results 

In this paper, we systematically study commuting Pauli Hamiltonians that are translation- 
invariant. We always assume that our Hamiltonians are frustration-free; every term in the 
Hamiltonian is minimized on the ground space. 

The main observation is that there is a purely algebraic description of commuting Pauli 
Hamiltonians in terms of maps between free modules over a Laurent polynomial ring by 
exploiting the translation-invariance. A Pauli matrix can be written as two binary numbers 
if we ignore the phase factors. For example, I = (00), cr^ = (10), = (01), = (11). 
We write these binary numbers in the coefficients of Laurent polynomials. The exponents 
of the Laurent polynomials will represent the positions at which the Pauli matrices act. 
If the Hamiltonian is translation-invariant, and there are finitely many distinct interaction 
types, then it follows that only a finite number of the Laurent polynomials convey all data 
of the Hamiltonian. We view this finite data as a map between two free modules over the 
translation-group algebra. We will show that the physical phase is solely determined by the 
image of this associated map. 

We provide a few tools to compute the transformations of the Hamiltonians by local 
unitary operators and coarse-graining. They come down to a well-defined set of elementary 
row operations on the matrices associated to the Hamiltonians. As we restrict our scope 
to the commuting Pauli Hamiltonians, the local unitary operators are also restricted to the 
Clifford operators. (Clifford operators maps a tensor product of Pauli operators to a tensor 
product of Pauli operators.) 

We define the characteristic dimension d associated to the Hamiltonian. If a Hamiltonian 
gives rise to a map between free modules, it is natural to think of the determinantal ideal of 
this map. The characteristic dimension is the KruU dimension of the algebraic set defined 
by this ideal. It is always upper bounded by the spatial dimension D. Moreover, d is less 
than or equal to Z) — 2 if the Hamiltonian is locally topologically ordered. 

The characteristic dimension d controls the rate at which the ground state degeneracy may 
increase. Roughly speaking, the logarithm of degeneracy can grow like L'^ where L is the 
linear system size. Thus, D = 3 is the minimal spatial dimension such that the degeneracy 
of a topologically ordered system can be diverging. For instance, the toric code models in 
various dimensions all correspond to the characteristic dimension 0, while the 3D cubic code 
model has characteristic dimension 1. However, it should be pointed out that the actual 
degeneracy does not behave as smooth as the function L'^; it can depend very sensitively 
on the system size. Indeed, it shall be shown that the degeneracy is related to the number 
of points in an algebraic set. The boundary condition imposes a global constraint on the 
relevancy of the points in the algebraic set. The numerically observed phenomena for the 
cubic code model shall be exactly calculated. 

We characterize topological charges in terms of torsion submodule. In particular, this 
characterization implies that a charge always appears as a vertex of a fractal operator, and 
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can be separated arbitrarily far from its partners by a local process with energy cost at most 
logarithmic in the separation distance. Here, the local process means a sequence of Pauli 
operators that are obtained by successive applications of single qubit Pauli operators. 

Using our formalism we derive a number of consequences of the low physical dimensionality. 
In one dimension, we algorithmically show how to transform an arbitrary Hamiltonian into 
several copies of the Ising models. In two dimensions, we characterize how the charges 
behave for topologically ordered models. Specifically, we prove that any excited state is a 
configuration of finitely many kinds of the charges, and the charges can be moved to an 
arbitrary position by a string operator. The result is a refined formulation of that of Bombin 
et al pm ITT] . In three dimensions, we prove that there always exists a point-like charge for 
any locally topologically ordered translation-invariant commuting Pauli Hamiltonian. This 
is a fundamental property of the three dimensions. It suggests that we might not be able 
to have a topologically ordered system in three dimensions where the excitations are all 
loop-like as in 4D toric code model. 

Our language is not completely new. A similar one appears in the error correcting code 
theory of computer science in the topic of multi-dimensional cyclic codes; see e.g. [12] and 
references therein. Also, there is an algebraic-geometry based design like Goppa codes 
However, the focus is different: We are interested in a fixed set of generators and exact 
sequences of modules describing the topological order. The fact that the commuting Pauli 
operators are represented as matrix, is very well-known in the theory of quantum error 
correcting codes. Our treatment is different in that the entries are not the binary values but 
the Laurent polynomials. 

Only the system of qubits, or spin-1/2, will be discussed, but all of our results and argu- 
ment straight-forwardly generalize to the system of qudits of prime dimensions. Technically, 
the ground field F2 for the qubit should be replaced by ¥p for a prime integer p. It is im- 
portant that the ground field is finite. Some numerical value 2 should be replace by the 
characteristic p of the field. With this generalization in mind, we keep necessary minus (— ) 
signs in the statements, which should be ignored for qu6its. Examples in this direction can 
be found in [T^ . 

We start by deriving the matrix representation of a commuting Pauli Hamiltonian, and 
explain in detail how the translation-invariance is exploited. The notion of modules over the 
translation-group algebra shall naturally emerge. Some operations on modules are induced 
by local unitary transformations on the physical system. They will define an equivalence 
relation between the matrix representations, and hence between Hamiltonians. We move on 
to the topological order and translate the conditions into those on a complex of modules. 
The characterization of charges will be given by exploiting the positive characteristic of the 
ground field. Finally, consequences of the topological order condition in two and three spatial 
dimensions will be derived. Explicit calculations and more examples are presented in the 
last section. All ring in the present paper shall be commutative with 1. 

Acknowledgements: The author would like to thank Sergey Bravyi, Lawrence Chung, 
Alexei Kitaev, John Preskill, Eric Rains, and Ari Turner for useful discussions. The author 
is supported in part by the Institute for Quantum Information and Matter, an NSF Physics 
Frontier Center, and the Korea Foundation for Advanced Studies. 
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binary field {0, 1} 
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2Fl, 
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number of qubits per site 

number of interaction types 

free i?-module of the interaction labels (rank t) 

free /^-module of Pauli operators (rank 2q) 

free i?-module of excitations (rank t) 

G ^ P, generating matrix or map for the stabilizer module 
P ^ E, generating matrix or map for excitations 
antipode map of the group algebra R. 
transpose followed by antipode map 

anti-symmetric 2q x 2q matrix f j 



Table 1. Reserved symbols. Any ring in this paper is commutative with 1. 



1. Algebraic structure of commuting Pauli Hamiltonians 
1.1. Pauli group as a vector space. The Pauh matrices 

J)' ''y^i^i 0')' '^^^(o -1) 

satisfy 

(TaCTb = iSabcC^c, {<^a, C^b} = ^Sab- 

Thus, the Pauli matrices together with scalars ±l,±i form a group under multiplication. 
Given a system of qubits, the set of all possible tensor products of the Pauli matrices form a 
group, where the group operation is the multiplication of operators. If the system is infinite, 
physically meaningful operators are those of finite support, i.e., acting on all but finitely 
many qubits by the identity. We shall only consider this Pauli group of finite support, and 
call it simply the Pauli group. An element of the Pauli group is called a Pauli operator. 

Since any two elements of the Pauli group either commute or anti-commute, ignoring the 
phase factor altogether, one obtains an abelian group. Moreover, since any element O of 
the Pauh group satisfies = ±7, An action of Z/2Z on Pauh group modulo phase factors 
P/{±1, ±i} is well-defined, by the rule n - O = where n E Z/2Z. For F2 = Z/2Z being a 
field, P/{±1, ±i} becomes a vector space over F2. The group of single qubit Pauli operators 
up to phase factors is identified with the two dimensional F2- vector space. If A is the index 
set of all qubits in the system, the whole Pauli group up to phase factors is the direct sum 
®ieA^ where Vi is the vector space of the Pauli operators for the qubit at i. Explicitly, 
/ = (00), (Tx = (10), = (01), = (11). A multi-qubit Pauli operator is written as a finite 
product of the single qubit Pauli operators, and hence is written as a binary string in which 
all but finitely many entries are zero. A pair of entries of the binary string describes a single 
qubit component in the tensor product expression. The multiplication of two Pauli operators 
corresponds to entry- wise addition of the two binary strings modulo 2. 
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The commutation relation may seem at first lost, but one can recover it by introducing a 

symplectic form [T^. Let -^i ~ ^ tie a symplectic form on the vector space (F2)^ of 

a single qubit Pauli operators. Q One can easily check that the commutation relation of two 
Pauli matrices Oi, O2 is precisely the value of this symplectic form evaluated on the pair of 
vectors representing Oi,02, respectively. Two multi-qubit Pauli operator (anti-) commutes 
if and only if there are (odd)even number of pairs of the anti- commuting single qubit Pauli 
operators in their tensor product expression. So, the two Pauli operator (anti-) commutes 
precisely when the value of the direct sum of symplectic form ©^g^ Ai is (non-)zero. (A 
could be infinite but the form is well-defined since any vector representing a Pauli operator 
is of finite support.) We shall call the value of the symplectic form the commutation value. 

1.2. Pauli space on a group. Let A be the index set of all qubits, and suppose now that 
A itself is an abelian group. There is a natural action of A on the Pauli group modulo phase 
factors induced from the group action of A on itself by multiplication. For example, if A = Z, 
the action of A is the translation on the one dimensional chain of qubits. li R = F2[A] is 
the group algebra with multiplicative identity denoted by 1, the Pauli group modulo phase 
factors acquires a structure of an i?-module. We shall call it the Pauli module. The Pauli 
module is free and has rank 2. 

Let r f be the antipode map of R, i.e., the F2-linear map into itself such that each 
group element is mapped to its inverse. Since A is abelian, the antipode map is an algebra- 
automorphism. Let the coefficient of a e i? at (7 G A be denoted by a^. Hence, a = X^^gA ^g9 
for any a E R. One may write ag = {ag)i. 

Define 

tr(a) = ai 

for any a E R. 

Proposition 1.1. Let {a,b),{c,d) G R^ be two vectors representing Pauli operators Oi,02 
up to phase factors: 




.geA / VgGA 

where a^^^ denotes the single qubit Pauli operator at g E A. Then, Oi and O2 commute if 
and only if 



Proof. The commutation value of (cri^^)"(cri^^)'^ and (o"i^'')"'(o"i^'')™'' is nm' — mn' E ¥2. 
Viewed as pairs of group algebra elements, ((Ti^^)"(o"i^'')™' and (cri^'')"'(cri^^)"*' are {ng,mg) 
and {n'g,m'g), respectively. We see that 



nm' — mn 



^The minus sign is not necessary for qubits, but is for qudits of prime dimensions 
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Since any Pauli operator is a finite product of tliese, tlie result follows by linearity. □ 

We wish to characterize a F2-subspace S of the Pauli module invariant under the action 
of A, i.e., a submodule, on which the commutation value is always zero. As we will see 
in the next subsection, this particular subspace yields a local Hamiltonian whose energy 
spectrum is exactly solvable, which is the main object of this paper. Let (a, b) be an element 
of 5* C i?^ = (F2[A])^. For any r & R, {ra,rb) must be a member of S. Demanding that the 
symplectic form on S vanish, by Proposition 11.11 we have 

tr(ra6 — rbd) = 0. 

Since r was arbitrary, we must have ab — ba = 0. Let us denote (a 6) as , and 

write any element of i?^ as a 2 x 1 matrix. We conclude that is a submodule of over 
R generated by si, . . . , such that any commutation value always vanishes, if and only if 

s\XiSj = 

for all i, j = 1, . . . ,t. 

The requirement that A be a group might be too restrictive. One may have a coarse 
group structure on A, the index set of all qubits. We consider the case that the index set is a 
product of a finite set and a group. By abuse of notation, we still write A to denote the group 
part, and insist that to each group element are associated q qubits (g > 1). Thus obtained 
Pauli module should now be identified with i?^"^, where R = F2[A] is the group algebra that 
encodes the notion of translation. We write an element v of R^'^ by a 2g x 1 matrix, and 
denote by the transpose matrix of v whose each entry is applied by the antipode map. 
We always order the entries of v such that the upper q entries describes the cXaj-part and the 
lower the cr^-part. Since the commutation value on R^'^ is the sum of commutation values 
on R^, we have the following: If is a submodule of i?^^ over R generated by Si, . . . , s^, the 
commutation value always vanishes on S, if and only if for alH, j = 1, . . . ,t 

sjXgSj = 

where = ^ ^^^^ is a 2g x 2q matrix. 

Let us summarize our discussion so far. 

Proposition 1.2. On a set of qubits Ax {1, . . . , g} where A is an abelian group, the group 
of all Pauli operators of finite support up to phase factors, form a free module P = R^'^ over 
the group algebra R = F2[A]. The commutation value 

(a, b) = tr(a^Aq6) 

for a,b & P is zero if and only if the Pauli operators corresponding to a and b commute. If a 
is a 2q X t matrix whose columns generate a submodule S ^ P, then the commutation value 
on S always vanishes if and only if 

a^Xqa = 0. 



'A symmetric bilinear form (r, s) = tr(rs) on R is non-degenerate. 
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1.3. Local Hamiltonians on groups. Recall that we place q qubits on each site of A. The 
total system of the qubits is A x {1, . . . , q}. 

Definition 1. Let 



be a local Hamiltonian consisted of Pauli operators that is (i) commuting, (ii) translation- 
invariant up to signs, and (iii) frustration-free. We call H a code Hamiltonian (also known 
as stabilizer Hamiltonian). The stabilizer module of H is the submodule of the Pauli module 
P generated by the images of /ii, . . . , /if in P. The number of interaction types is t. 

The energy spectrum of the code Hamiltonian is trivial; it is discrete and equally spaced. 

Example 1. One dimensional Ising model is the Hamiltonian 



igZ 

The lattice is the additive group Z, and the group algebra is R = F2[a;, x]. The Pauli module 
is and the stabilizer module S is generated by 



1.4. Excitations. For a code Hamiltonian if, an excited state is described by the terms 
in the Hamiltonian that have eigenvalues —1. Each of the flipped terms is interpreted as 
an excitation. Although the actual set of all possible configurations of excitations that are 
obtained by applying some operator to a ground state, may be quite restricted, it shall be 
convenient to think of a larger set. Let E be the set of all configurations of finite number of 
excitations without asking physical relevance. Since an excitation is by definition a flipped 
term in H, the set E is equal to the collection of all finite sets consisted of the terms in H. 

If Pauh operators Ui,U2 acting on a ground state creates excitations 61,62 G E, their 
product U1U2 creates excitations (61 U 62) \ (61 fl 62). Here, we had to remove the intersection 
because each excitation is its own annihilator; any term in the H squares to the identity. 
Exploiting this fact, we make E into a vector space over F2. Namely, we take formal linear 
combinations of terms in H with the coefficient 1 G F2 when the terms has —1 eigenvalue, 
and the coefficient G F2 when the term has -|-1 eigenvalue. The symmetric difference is now 
expressed as the sum of two vectors 61 + 62 over F2. In view of Pauli group as a vector space, 
U1U2 is the sum of the two vectors vi + V2 that respectively represents Ui,U2. Therefore, 
the association f/j 1— )■ 6j induces a linear map from the Pauli space to the space of virtual 
excitations E. 

The set of all excited states obeys the translation-invariance as the code Hamiltonian 
H does. So, is a module over the group algebra R = F2[A]. The association Ui 6j 
clearly respects this translation structure. Our discussion is summarized by saying that the 
excitations are described by an i?-linear map 






One can view this as the matrix a of Proposition 11.21 H is commuting; cr^Xicr = 0. 







e:P^E 



from the Pauli module P to the module of virtual excitations E. 
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As the excitation module is the collection of all finite sets of the terms in H , we can speak 
of the module of generator labels G, which is equal to E as an i?-module. G is a free module 
of rank t if there are t types of interaction. The matrix a introduced in Section 11.21 can be 
viewed as 

from the module of generator labels to the Pauli module. 

Proposition 1.3. If a is the generating map for the stabilizer module of a code Hamiltonian, 
then 

e = a^Xg. 

Proof. This is a simple corollary of Proposition 11.21 Let hi^g be the terms in the Hamiltonian 
where i = 1, . . . ,t, and g ^ A. In the Pauli module, they are expressed as ghi where hi is 
the 2-th column of cr. For any u & P, let e{u)i be the i-th component of ^{u). By definition, 

= X^S- {ighi)^Ku) = X^S- (ghlXguj = h\\qU 

geA geA 

Thus, hlXg is the i-th row of e. □ 

Remark 1. The commutativity condition in Proposition 11.21 of the code Hamiltonian is 
recast into the condition that 

G^P^E 
be a complex, i.e., e o a = 0. Equivalently, 

im 0" C (im o") ^ = ker e 
where _L is with respect to the symplectic form. 

2. Equivalent Hamiltonians 

The stabilizer module entirely determines the physical phase of the code Hamiltonian in 
the following sense. 

Proposition 2.1. Let H and H' be code Hamiltonians on a system of qubits, and suppose 
their stabilizer modules are the same. Then, there exists a unitary 

geA 

mapping the ground space of H onto that of H' . Moreover, there exist a continuous one- 
parameter family of gapped Hamiltonians connecting UHW and H' . 

Proof. Let {pa} be a maximal set of F2-linearly independent Pauli operators of finite support 
that generates the common stabilizer module S. {pa} is not necessarily translation-invariant. 
Any ground state of if is a common eigenspace of {pa} with eigenvalues p^ \ip) = I?/'), 
Ca = ±1. Similarly, the ground space of H' gives the eigenvalues e'„ = ±1 for each p^. 

The abelian group generated by {pa} is precisely the vector space S*, and the assignment 
Pa ^ Ca defines a dual vector on S. If f/ is a Pauli operator of possibly infinite support, 
then Pali \ip) = e'^CaU \ip) for some e'^ = ±1, where e'^ is determined by the commutation 
relation between U and pa- Thus, the first statement follows if we can find U such that 
the commutation value between U and pa is precisely e'^. This is always possible since the 
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dual space of the vector space P is isomorphic to the direct product HaxIi <j}^25 which is 
vector-space-isomorphic to the Pauh group of arbitrary support up to phase factors. H 

Now, UHW and H' have the same eigenspaces, and in particular, the same ground space. 
Consider a continuous family of Hamiltonians 

H{u, u) = uUHU^ + uH' 

where u, u' G M. It is clear that 

H = H{1, 0) ^ H{1, 1) ^ H{0, 1) = H' 
is a desired path. □ 

The criterion of Proposition 12.11 to classify the physical phases is too narrow. Physically 
meaningful universal properties should be invariant under simple and local changes of the 
system. More concretely. 

Definition 2. Two code Hamiltonians H and H' are equivalent if their stabilizer modules 
become the same under a finite composition of symplectic transformations, coarse-graining, 
and tensoring ancillas. 

We shall define the symplectic transformations, the coarse-graining, and the tensoring 
ancillas shortly. 

2.1. Symplectic transformations. 

Definition 3. A symplectic transformation T is an automorphism of the Pauli module 
induced by a unitary operator on the system of qubits such that 

TU,T = A, 

where f is the transposition followed by the entry-wise antipode map. 

Only the unitary operator on the physical Hilbert space that respects the translation can 
induce a symplectic transformation. By definition, a symplectic transformation maps each 
local Pauli operator to a local Pauli operator, and preserves the commutation value for any 
pair of Pauli operators. 

Proposition 2.2. Any two unitary operators Ui, U2 that induces the same symplectic trans- 
formation differ by a Pauli operator (of possibly infinite support). 

Proof. The symplectic transformation induced by U = UIU2 is the identity. Hence, U 
maps each single qubit Pauli operator ai^^z'' to ±ai?j\ By the argument as in the proof of 
Proposition 12. H there exists a Pauli operator O of possibly infinite support that acts the 
same as U on the system of qubits. Since Pauli operators form a basis of the operator algebra 
of qubits, we have O = U. □ 

The effect of a symplectic transformation on the generating map a is a matrix multiplica- 
tion on the left. 

a Ua 



If F is a finite dimensional vector space over some field, the dual vector space of 0j V is isomorphic to 
y where I is an arbitrary index set. 
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For example, the following is induced by uniform Hadamard, controUed-Phase, and controUed- 
NOT gates. For notational clarity, define -E'ij(a) (i ^ j) as the row-addition elementary 
2q X 2q matrix 

where S^,, is the Kronecker delta and a E R — ¥2[A]. Recall that we order the components 
of P such that the first half components are for cr-r-part, and the second half components are 
for (Jz-part. 

Definition 4. The following are elementary symplectic transformations: 

• (Hadamard) Ei^i+g{-l)Ei+g^i{l)Ei^i^g{-l) where 1 <i <q, 

• (controUed-Phase) Ei^q^i{f ) where f = f and I < i < q, 

• (controUed-NOT) Eij{a)Ej+q^i+g{—d) where I < i ^ j < q. 

Recall that the Hadamard gate is a unitary transformation on a qubit given by 

with respect to basis {|0) , |1)}. At operator level, 

UhXUI = Z, UhZUI, = X 

where X and Z are the Pauli matrices and a^, respectively. Thus, the application of 
Hadamard gate on every i-th qubit of each site of A swaps the corresponding X and Z 
components of P. 

The controlled phase gate is a two-qubit unitary operator whose matrix is 



/I 








M 





1 














1 





^0 








-1/ 



with respect to basis {|00) , |01) , |10) , |11)}. At operator level, 

Up(X ® I)Ul = X ® Z, Up{Z®I)Ul = Z®I, 

Up{I ®x)ul = Z ®X, Up{I ® z)ul = I ® z. 

Note that since f/p is diagonal, any two Up on different pairs of qubits commute. Let [g-ii) 
denote the i-th qubit at gf G A. The uniform application 

U^'' = X{Up{{Ki),{h + g,i)) 

heK 

of Up throughout the lattice A such that each Up{{h, i), {h + g, i)) acts on the pair of qubits 
{h,i) and {h + g,i) is well-defined. From the operator level calculation of Up, we see that 
Ug^ induces 

P 3 {. . . ,Xi, . . . ,Zi, . . .) {. . . ,Xi, . . . ,Zi + {g + g)xi, . . .) e P 
on the Pauli module, which is represented as Ei^gi{g + g). The composition 

Tri^TTi^ . . . Tr(.^ 

gi 92 ^ Qn 
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of finitely many controlled-Phase gates Ug with different g is represented as Ei^g^i{f) where 
f = f = Ylk=i 9k+gk- The single qubit phase gate 




maps X i-^ Y and Z ^ Z. On the Pauli module P, it is 

-P 9 (. . . , Xj, . . . , 2^j, . . .) I ^ (. . . , Xj, . . . , -|- Xj, . . .) (z p. 

which is Ei+g^i{l). Note that any f E R such that / = / is always of form f = J2 9k + Ok or 
f = 1 + Yl 9k + dk where Qk are monomials. Thus, the Phase gate and the controlled-Phase 
gate induce transformations Ei+g^i{f) where f = f- 

The controlled-NOT gate is a two-qubit unitary operator whose matrix is 



/I 








o\ 





1 

















1 


^0 





1 





with respect to basis {|00) , |01) , |10) , |11)}. That is, it flips the target qubit conditioned on 
the control qubit. At operator level, 

Un{X^I)uI = X^X, Un{Z®I)uIj = Z®I, 

Un{i®x)uI = i®x, UN{I®Z)Ulf = z®z. 

\i i < j, the uniform application 

U^'-^^ = (^Up{{h,z),{h + g,j)) 

such that each UN{{h,i), [h + g,])) acts on the pair of qubits {h^i) and {h + g^j) with one 
at (/i, i) being the control induces 

P 9 ( . . . , X j , . . . , Xj , . . . , , . . . , , . . . ) 

I y (. . . , Xj, . . . , Xj -|- gxi^ . . . , Zi -\- gzj, . . . , zj , . . G P. 

Thus, any flnite composition of controlled-NOT gates with various g is of form Eij{a)Ej^q i^g{d). 
It might be useful to note that the controlled-NOT and the Hadamard combined, induces a 
symplectic transformation 

• (controlled-NOT-Hadamard) Ei^qj{a)Ej^q^i{d) where a & R and 1 < i ^ j < q. 
Remark that an arbitrary row operation on the upper q components can be compensated by 
a suitable row operation on the lower q components so as to be a symplectic transformation. 

2.2. Coarse-graining. Not all unitary operators conform with the lattice translation. In 
Example [1] the lattice translation has period 1. Then, for example, the Hadamard gate 
on every second qubit does not respect this translation structure; it only respects a coarse 
version of the original translation. We need to shrink the translation group to treat such a 
unitary operators. 

Let A be the original translation group of the lattice with q qubits per site, and A' be 
its subgroup of flnite index: |A/A'| = c < oo. The total set of qubits A x {1, . . . , g} is 
set-theoretically the same as A' x {1, . . . , c} x {1, . . . , g} = A' x {1, . . . , cq}. We take A' 
as our new translation group under coarse-graining. The Pauli group modulo phase factors 
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remains the same as a F2-vector space for it depends only on the total index set of qubits. 
We shall say that the system is coarse-grained by R' = F2[A'] if we restrict the scalar ring R 
to R' for all modules pertaining to the system. 

For example, suppose A = Z^, so the original base ring is R = ¥2[x,y,x,y]. If we coarse- 
grain by R' = ¥2[x', y',x', y'] where x' = x"^, y' 
original lattice as a single new site. 



y^, we are taking the sites 1, x, ?/, xy of the 



2.3. Tensoring ancillas. We have considered possible transformations on the stabilizer 
modules of code Hamiltonians, and kept the underlying index set of qubits invariant. It is 
quite natural to allow tensoring ancilla qubits in trivial states. In terms of the stabilizer 
module 5* C P = i?^'?, it amounts to embed S into the larger module i?^'' where q' > q. 

' be the generating matrix of S as in Proposition 11.21 By tensoring 



Concretely, let a 
ancilla, we embed S as 



crx 





o\ 








o-z 







V 



This amounts to taking the direct sum of the original complex 

G^P^E 

and the trivial complex 



O^R 



^R' 



1 0) 



> i?^ 



to form 



G®R^P®R^^E®R. 



3. Topological order 

From now on we assume that A is isomorphic to as an additive group. D shall be 
called the spatial dimension of A. 

Definition 5. Let a : G — ?► P be the generating map for the stabilizer module of a code 
Hamiltonian H. We say H is exact if (ima)^ = ima, or equivalently 

is exact, i.e., kere = imo". 

It follows that the exactness condition is a property of the equivalence class of code Hamil- 
tonians in the sense of Definition [2j 

By imposing periodic boundary conditions, a translation-invariant Hamiltonian yields a 
family of Hamiltonians {H{L)} defined on a finite system consisted of sites. One might be 
concerned that some H{L) would be frustrated. We intentionally exclude such a situation. 
The frustration might indeed occur, but it can easily be resolved by choosing the signs 
of terms in the Hamiltonian. In this way, one might loose the translation-invariance in a 
strict sense. However, we retain the physical phase regardless of the sign choice because 
different sign choices are related by a Pauli operator acting on the whole system which is a 
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product unitary operator. Hence, the entanglement property of the ground state and the all 
properties of excitations do not change. 

Definition 6. Let H{L) be Hamiltonians on a finite system of linear size L in D dimen- 
sional physical space, and Ul be the corresponding ground space projector. H{L) is called 
topologically ordered if for any O supported inside a hypercube of size one has 

(1) HlOUl oc n^. 

This means that no local operator is capable of distinguishing different ground states. 
This condition is trivially satisfied if H{L) has a unique ground state. A technical condition 
that is used in the proof of the stability of topological order against small perturbations is 
the following 'local topological order' condition p^HTlIIB]- We use a simplified version. We 
say pyramid region A{r) of linear size r for the set 




or its translation. The apex will be referred to the point of the smallest coordinates. 

Definition 7. Let H{L) be code Hamiltonians on a finite system of linear size L in D 
dimensional physical space. For any pyramid region A = A{r) of linear size r, let H^ be the 
projector onto the common eigenspace of the smallest eigenvalue of terms in the Hamiltonian 
H{L) that are supported on A. For 6 > 0, denote by A^ the distance b neighborhood of A. 
H{L) is called locally topologically ordered if there exists a constant 6 > such that for any 
operator O supported on a pyramid region A of linear size r < L/2 one has 

(2) U^bOUAb oc H^i,. 

Since any operator is a C-linear combination of Pauli operators, if Eq. dl]),® are satisfied 
for Pauli operators, then (local) topological order condition follows. If a Pauli operator O is 
ant i- commuting with a term in a code Hamiltonian H{L), The left-hand side of Eq. ([I]),© are 
identically zero. In this case, there is nothing to be checked. If O acting on A is commuting 
with every term in H{L) supported inside A^, Eq. ([1]) demands that it act as identity on 
the ground space, i.e., O must be a product of terms in H{L) up to ±i, ±1. Eq. ([2]) further 
demands that O must be a product of terms in H{L) supported inside A'' up to ±i, ±1. 

Lemma 3.1. A code Hamiltonian H is exact if and only if H{L) is locally topologically 
ordered for all sufficiently large L. 

It shall be important to use Laurent polynomials to express elements of the group algebra 
R = ¥2[Z^] = F2[xi, x]""'", • • • , xd, x]^^]. For example, 

xy^z'^+xy^^ ^ 1(1,2,2) + 1(1,-1,0). 

The sum of the absolute values of exponents of a monomial will be referred to as absolute 
degree. The absolute degree of a Laurent polynomial is defined to be the maximum absolute 
degree of its terms. The degree measures the distance or size in the lattice. 

Proof. ( "only if" ) Without loss of generality, assume that a is expressed in polynomials with 
non-negative exponents. Let w be the maximum of all absolute degree of entries of a; it is 
the interaction range. Let O be a Pauli operator supported on a pyramid region of linear 
size r. Choose the origin of the lattice such that the pyramid region has its apex at 
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{w,w, . . . ,w) G . Then O as an element of the Pauh module is expressed as a 2g x 1 
matrix v whose entries are polynomials in D variables with positive exponents of absolute 
degree at most r + w. 

We have to show that for any such v, if v G kere, then v can be expressed as a linear 
combination 



of the columns at of a such that the coefficients ci G R have absolute degree not exceeding 
w + r. (The constant b in the definition of the local topological order will be 0{w).) 

By the exactness assumption, kere = imcr. Since all of our Laurent polynomials have non- 
negative exponents, we can use a tool for modules over polynomial rings: Grobner basis fi9\ 
Chapter 15]. Buchberger's algorithm to calculate Grobner basis of imcr with respect to the 
homogeneous degree monomial order produces Grobner basis of degree not exceeding w. The 
standard division algorithm can then be used to compute the coefficients q. Thus obtained 
Ci have non-negative exponents that are bounded above hj r + w. 

("if") Suppose V G kere. We have to show v G imcr. Choose so large L that the Pauli 
operator O representing v is contained in a pyramid region far from the boundary. The local 
topological order condition implies that O is a product of terms near the pyramid region. 
Since this product expression is independent of the boundary, we see f G im cr. □ 

The quoted theorem below characterizes an exact sequence from the properties of connect- 
ing maps. A few notions should be recalled. Let M be a matrix, not necessarily square, over 
a ring. A minor is the determinant of square submatrix of M. k-th determinantal ideal /a:(M) 
is the ideal generated by all A; x A; minors of M. It is not hard to see that the determinantal 
ideal is invariant under any invertible matrix multiplication on either side. The rank of M is 
the largest k such that fc-th determinantal ideal is nonzero. (Sometimes 0-th determinantal 
ideal is taken to be the unit ideal by convention.) For a map between free modules, we 
write J(0) to denote the fc-th determinantal ideal of the matrix of (j) where k is the rank of 
that matrix. Fitting Lemma [191 Corollary-Definition 20.4] states that determinantal ideals 
only depend on coker 0. 

The (Krull) dimension of a ring is the supremum of lengths of chains of prime ideals. 
Here, the length of a chain of prime ideals 



is defined to be n. In particular, in a dimension zero ring every prime ideal is maximal. 
The codimension or height of a prime ideal p is the supremum of the lengths of chains of 
prime ideals contained in p. That is, the codimension of p is the Krull dimension of the local 
ring Rp. The codimension of an arbitrary ideal / is the minimum of codimensions of primes 
that contain /. The dimension of an ideal / C i? is the Krull dimension of R/I. If S is 
an affine domain, i.e., a homomorphic image of a polynomial ring over a field with finitely 

^This part can be adapted to an error correcting procedure or a decoder. The bottleneck of the universal 
decoder presented in [9] is the routine that tests whether a given cluster of excitations can be created by a 
Pauli operator supported in the box that envelops the cluster. The Grobner basis for im e in the homogeneous 
degree monomial order provides a fast algorithm for it: The division algorithm yields zero remainder with 
respect to the Grobner basis, if and only if the given cluster is in im e. One minor modification is that one 
has to consider an enveloping pyramid in place of the enveloping box. Note also that this argument proves 
that the topological order condition as defined in is always satisfied if the code Hamiltonian is exact. 




Po C pi C •■■ C p 
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many variables such that S has no zero-divisors, it holds that codim / + dim / = dim 5* [T9t 
Chapter 13]. 

We shall be dealing with three different kinds of 'dimensions': The first one is the spatial 
dimension D, which has an obvious physical meaning. The second one is the Krull dimension 
of a ring, just introduced. The Krull dimension is upper bounded by the spatial dimension 
in any case. The last one is the dimension of some module as a vector space. Recall that all 
of our base ring contains a field - F2 for qubits. The vector space dimension arises naturally 
when we actually count the number of orthogonal ground states. The dimension as a vector 
space will always be denoted with a subscript like dim^j. 

Proposition 3.2. |19l Theorem 20.9, Proposition 18.2], [201 Chapter 6 Theorem 15] // a 

complex of free modules over a ring 

O^F„^F„_i^---^Fi^Fo 

is exact, then 

• rank = rank + rank for k = 1, . . . ,n — 1 

• rank F„ = rank 0„ . 

• = (1) or else codim 7(0^) > k for k = 1, . . . ,n. 

Remark 2. For exact code Hamiltonian, we have a exact sequence G ^ P E. 
As we will see in Lemma 17. H coker a has a finite free resolution, and we may apply the 
Proposition 13.21 Since /^(cr) = /^(e) for any A; > 1, we have 

2q = rank P = rank cr + rank e = 2 rank a. 

The size 2g x t of the matrix cr satisfies t > q. If Iq{cr) 7^ R, then codim /^((t) > 2. 

4. Ground state degeneracy 

Let H{L) be the Hamiltonians on finite systems obtained by imposing periodic boundary 
conditions as in Section [31 A symmetry operator of H{L) is a linear combination of Pauli 
operator that commutes with H{L). In order for a Pauli symmetry operator to have a 
nontrivial action on the ground space, it must not be a product of terms in H{L). In 
addition, since H{L) is a sum of Pauli operators, a symmetry Pauli operator must commute 
with each term in H{L). Hence, a symmetry Pauli operator O with nontrivial action on the 
ground space must have image v in the Pauli module such that 

v{0) G ker \ im ctx 

where 

G/bLG ^ P/biF ^ E/bLE 

and 

bi = (a;f-l,...,4-l)Ci?, 

which effectively imposes the periodic boundary conditions. Since each term in H{E) acts 
as an identity on the ground space, if O' is a term in H{L), the symmetry operator O and 
the product 00' has the same action on the ground space. 00' is expressed in the Pauli 
module as v{0) + v'{0') for some v' G imai- Therefore, the set of Pauli operators of distinct 
actions on the ground space is in one-to-one correspondence with the factor module 

K{L) =kereL / imai- 
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The vector space dimension diniFa K{L) is precisely the number of independent Pauh oper- 
ators that have nontrivial action on the ground space. Since kere^, = (imcr^,)-'- by definition 
of e, and imcr^ as an F2-vector space is a null space of the symplectic vector space P/b^P, 
it follows that {im aj^)-^ = imai^W for some hyperbolic subspace W. The quotient space 
K{L) = W is thus hyperbolic and has even vector space dimension 2k. Choosing a symplec- 
tic basis for K{L), it is clear that K{L) represents the tensor product of k qubit- algebras. 
Therefore, the ground space degeneracy is exactly 2'^. In the theory of quantum error cor- 
recting codes, k is called the number of logical qubits, and the elements of K{L) are called 
the logical operators. In this section, k will always denote | dim^^ K. 

4.1. Condition for degenerate Hamiltonians. 

Definition 8. The associated ideal for a code Hamiltonian is the g-th determinantal ideal 
Iq{(^) C P of the generating map a. Here, q is the number of qubits per site. The charac- 
teristic dimension is the KruU dimension of R/Iq{a). 

Lemma 4.1. Let I he the associated ideal of an exact code Hamiltonian, and mhe a maximal 
ideal of P. Then, / ^ m implies that the localized homology 



is zero for all L > 1. 

Proof. Recall that the localization and the factoring commute. By assumption, (/g(e))m = 

{Iq{a))m = (1) = Pm =: S. For notational convenience in the proof, let us drop the subscript 
m denoting the localization. Recall that the local ring 5* has the unique maximal ideal 
m, and any element outside the maximal ideal is a unit. If every entry of e is in m, then 
Iq{^) ^ m 7^ S". Therefore, there is a unit entry, and by column and row operations, e is 
brought to 



where e' is the submatrix. It is clear that /g_i(e') C Iq{e) since any g — 1 x g — 1 submatrix of 
e' can be thought of as a g x g submatrix of e where the first column and first row have the 
unique nonzero entry 1 at (1, 1). It is also clear that /q_i(e') ^ Ig{e) since any g x g submatrix 
of e contains either zero row or column, or the (1, 1) entry 1 of e. Hence, /q_i(e') = (1), and 
we can keep extracting unit elements into the diagonal by row and column operations ^201 
Chapter 1 Theorem 12]. After g steps, t x 2q matrix e becomes precisely 



where idg is the q x q identity matrix. Since localization preserves the exact sequence 
G ^ P ^ E, a maps to the lower g components of P with respect to the basis where e is in 
the above form. Since Iq{cr) = (1), we must have (after basis change) 



Therefore, even after factoring by the proper ideal bi, the homology K{L) = kere^ / imai 
is still zero. □ 



P'(L)m = ker(eL)n, / im(o-L) 



m 
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Corollary 4.2. The associated ideal of an exact code Hamiltonian is the unit ideal, i.e., 
Iq{(y) = R, if and only if 

K{L) = kere^ / imcTL = 

for all L > 1. 

Proof. If I{a) = R, I{a) is not contained in any maximal ideal m. The above lemma says 
K{L)^-n = 0. Since a module is zero if and only if its localization at every maximal ideal is 
zero, K{L) = for all L > 1. 

For the converse, observe that if F is any extension field of F2, for any F2-vector space W, 
we have dim^ F W = dim^^ W. We replace the ground field F2 with its algebraic closure 
F" to test whether K{L) 7^ 0. If Iq{cr) is not the unit ideal, then it is contained in a maximal 
ideal m C. R. By Nullstellensatz, m = {xi — ai, . . . ,xd — cld) for some G F". Since in R 
any monomial is a unit, we have 7^ 0. Therefore, there exists L > 1 such that af = 1 and 
2 \ L. The equation x^ — 1 = has no multiple root. 

We claim that K{L) 7^ 0. It is enough to verify this for the localization at m. Since 
anything outside m is a unit in R^ and each xf — 1 contains exactly one Xj — Oj factor, we 
see {bL)m = vrim- Therefore, {eL)m = em/(^'L)m and ((Tl)^ = cr^/{bL)m is a matrix over the 
field R/m = F". Since Ig{<j) C m, we have Ig{crL)m = 0. That is, rank]Fa(cr£^)^ < q and 
Tajikfa(eL)m < q- It is clear that dimpa K^L)^ = dim^a ker(e2,)n,/ im((Ti)m > 2. □ 

This corollary says that in order to have a degenerate Hamiltonian H{L), one must have 
a proper associated ideal. We shall simply speak of a degenerate code Hamiltonian if its 
associated ideal is proper. 

4.2. Counting points in an algebraic set. It is important that the factor ring 

R/bL = F2[xi, . . . , x^] / (xf - 1, . . . , x£ - 1) 

is finite dimensional as a vector space over F2, and hence is Artinian. In fact, dimpj R/bi = 
L^. Due to the following structure theorem of Artinian rings, K{L) can be explicitly analyzed 
by the localizations. 

Proposition 4.3. [21, Chapter 8] [19, Section 2.4] Let S be an Artinian ring. (For example, 
S is a homomorphic image of a polynomial ring over finitely many variables with coefficients 
in a field ¥, and is finite dimensional as a vector space over ¥.) Then, there are only finitely 
many maximal ideals of S , and 

m 

where the sum is over all maximal ideals xn of S and is the localization of S at m. 

The following calculation tool is sometimes useful. Recall that a group algebra is equipped 
with a non-degenerate scalar product {v, w) = tr^vw). This scalar product naturally extends 
to a direct sum of group algebras. 

Lemma 4.4. Let ¥ be a field, and S = F[A] be the group algebra of a finite abelian group A. 
Let Vi be elements of a free S-module S", and 




= {ve S'''\{v,v^) = for alii}. 
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Then, the dual vector space N* of N is isomorphic as a vector space to 

N* = S''/J2Svi. 

i 

Proof. Consider : S*" 9 x i— )■ (•, x) G N*. The map (f) is surjective since the scalar product 
is non-degenerate and is a finite dimensional vector space. The kernel of is precisely 
N^. □ 

Corollary 4.5. Put 2k = dim^^ K{L). Then, 

k = qL^ — dimpj im ai = dimp^ ker — qL^. 

Further, if q = t, then 

k = dimpj coker ei,. 



Proof. Put 5* = R/hi- If f i, . . . , Vt denote the columns of ctl, we have 
(3) ker a\ = Xg ker = P| = Sv-^ = (im ctl)"^ • 

Hence, dimpa ker ez, = dimpj S'^''— dimpj imai. Since dimpa S = and K{L) = kere^/ im ctl, 
the first claim follows. 

Since imaL = 5"*/ ker ctl, if t = q, we have k = dimpj ker ctl by the first claim. From 
Eq. ([3]), we conclude that k = dimpj 5**/ imcrj^ = dimpj cokere^. □ 

We will apply these results in Section O 

The characteristic dimension is related to the rate at which the degeneracy increases as 
the system size increases in the following sense. Recall that 2k = dimpj K{L) and the ground 
state degeneracy is 2^^. 

Lemma 4.6. Suppose 2 \ L. Let he the algebraic closure of ¥2. If N is the number of 
maximal ideals in ^^^^ contains bi + Iq{cr), then 

N < dimp^ K{L) < 2qN. 



Proof. We replace the ground field F2 with F". Any maximal ideal of an Artinian ring 
F"[x^^]/bi is of form m = (a;i — ai, . . . , — od) where af = 1 by Nullstellensatz. Since 
2 I L, we see that {bL)m = xn^ and that {R/bL)m = F'* is the ground field. (See the proof of 
Corollary 1121) 

Now, Iq{cr) + bi C m iff Iq{a)m + {bL)m ^ = (^'L)m iff Iq{o') becomes zero over 
-Rm/(&L)m = F" iff 1 < dimpa K{L)m < 2q. Since by Proposition 14. 3[ dimpa K{L) is a finite 
direct sum of localized ones, we are done. □ 

Lemma 4.7. Let I be an ideal such that dimR/I = d. We have 

dimF2 R/{I + bi) < cL^ 



for all L > 1 and some constant c independent of L. 
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Proof. We replace the ground field with its algebraic closure F". Write Xi for the image of Xi 
in R/I. By Noether normalization theorem Theorem 13.3], there exist yi, ■ ■ ■ ,yd ^ R/I 
such that R/I is a finitely generated module over F"[?/i, . . . ,y^. Moreover, one can choose 
Hi ~ XljLi ^ij^j for some rank d matrix M whose entries are in F". Making M into the 
reduced row echelon form, we may assume yi = Xi + J2j>d '^ij^j ^ach 1 <i < d. 

Let S = W^lzi, . . . , zd] be a polynomial ring in D variables. Let : 5* —)■ /?/(/+ bi,) be 
the ring homomorphism such that Zi ^ yi for 1 < i < d and Zj H- Xj for d < j < D. By the 
choice of yi, is clearly surjective. Consider the ideal J of S" generated by the initial terms 
of ker0 with respect to the lexicographical monomial order in which zi < ■ ■ ■ < Z£,. Since 
Xj is integral over F[?/i, . . . , t/^], the monomial ideal J contains z^^ for some positive rij for 
d < j < D. Here, nj is independent of L. Since G J for 1 < z < d, we conclude that 

dim^a R/ (J + bi) = dimpa S/ J < L'^ ■ 71^+112^+2 ■ ■ ■ 't-d 
by Macaulay theorem [19, Theorem 15.3]. □ 

Corollary 4.8. If 2 \ L, and d = dimi?/Jg(cr) is the characteristic dimension of a code 
Hamiltonian, then 

dimf, K{L) < cL'^ 
for some constant c independent of L. 

Proof. If J = bi + /(o"), in Lemma 14.61 is equal to dimpa R/ rad J. This is at most 
diiRfa R/ J. □ 

Lemma 4.9. Let d be the characteristic dimension. There exists an infinite set of integers 
{Li} such that 

dimp, K{Li) > Li72 

Proof. We replace the ground field with its algebraic closure F'*. Let p' 3 I {a) be a prime 
of R of codimension D — d. Let p be the contraction (pull-back) of p' in the polynomial ring 
S = ¥"-[xi, . . . ,xd]- Since the set of all primes of R is in one-to-one correspondence with 
the set of primes in S that does not include monomials, it follows that p has codimension 
D — d and does not contain any monomials. Let V denote the affine variety defined by 
p = {gi, . . . , gn). Since p contains no monomials, V is not contained in any hyperplanes 
Xi = {t = l,...,D). 

Let Ai be a finite subfield of F" that contains all the coefficients of gi, so V can be defined 
over Ai. Let An C F" be the finite extension fields of Ai of extension degree n. Put 
Ln = \An\ — 1. For any subfield A of F°, let us say a point of V is rational over A if its 
coordinates are in A. The number N'{Ln) of points (oj) G V satisfying af" = 1 is precisely 
the number of the rational points of V over An that are not contained in the hyperplanes 
Xi = 0. Since I (a) C p', the number N in Lemma [4.61 is at least N'{Ln)- It remains to show 
2iV'(L„) > for all sufficiently large n. 

This follows from the result by Lang and Weil [22], which states that the number of points 
of a projective variety of dimension d that are rational over a finite field of m elements is 

+ O ^m'^~^ j asymptotically in m. Since Lang- Weil theorem is for projective variety and 

we are with an affine variety V, we need to subtract the number of points in the hyperplanes 
Xi = {i = 0,1, . . . , D) from the Zariski closure of V. The subvarieties in the hyperplanes, 
being closed, have strictly smaller dimensions, and we are done. □ 
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5. Fractal operators and topological charges 



This section is to provide a characterization of topological charges, and their dynamical 
properties. Before we turn to general characterization and define fractal operators, let us 
review familiar examples. Note that for two dimensions the base ring is R = ¥2[x,x,y,y]. 

Example 2 (Toric Code). Although the original two-dimensional toric code has qubits on 
edges P , we put two qubits per site of the square lattice to fit it into our setting. Concretely, 
the first qubit to each site represents the one on its east edge, and the second qubit the one 
on its north edge. With this convention, the Hamiltonian is the negative sum of the following 
two types of interactions: 



XI XX 

I I 

II -IX 



ZI-II 

I I 

ZZIZ 



y-xy 
I I 
1 — X 



where we used X, Z to abbreviate ax,<Jz, and omitted the tensor product symbol. Here, the 
third square specifies the coordinate system of the square lattice. Since there are q = 2 qubits 
per site, the Pauli module is of rank 4. The corresponding generating map a : i?^ — )■ i?^ is 
given by the matrix 

f y + xy \ fl + x 

x + xy ^ l + y 



C"2D- toric 



v 













V 








l+y 
l + xj 

Here, the each column expresses each type of interaction. It is clear that 



l+y 1 + X 



e2D-toric — ^"^^2 



1 + X 





l + y 
1 + xJ 



l + y 




and kere = ima; the two dimensional toric code satisfies our exactness condition. The 
associated ideal is /(a) = ((1 + x)^, (1 + x)(l + y), (1 + y^)- The characteristic dimension 
is dim R/I{a) = 0. Note also that anncokere = {x — l,y — 1). The electric and magnetic 

charge are represented by G-E\ime, respectively. 



OJ ' ^1^ 

The connection with cellular homology should be mentioned, a can be viewed as the 
boundary map from the free module of all 2-cells with Z2 coefficients of the cell structure 
of 2-torus induced from the tessellation by the square lattice. Then, e is interpreted as the 
boundary map from the free module of all 1-cells to that of all 0-cells. cr or e is actually 
the direct sum of two boundary maps. Indeed, the space K{L) = kere/,/ imaz, of operators 
acting on the ground space (logical operators) has four generators 
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\l + y + ■■■ + 2/^-7 
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which correspond to the usual nontrivial first homology classes of 2-torus. 

The description by the cellular homology might be advantageous for the toric code over 
our description with pure Laurent polynomials; in this way, it is clear that the toric code 
can be defined on an arbitrary tessellation of compact orientable surfaces. However, it is 
unclear whether this cellular homology description is possible after all for other topologically 
ordered code Hamiltonians. 

Example 3 (2D Ising model on square lattice). The Ising model has nearest neighbor 
interactions that are horizontal and vertical. In our formalism, they are represented as 1 + x 
and l + y. Thus, 


1 + X 1 + y. 

As it is not topologically ordered, the complex G — > P — > £^ is not exact. Moreover, a is not 
injective. 

'l + y' 



C" 2D Ising 



C2D Isingjl 



1 + X 



generates the kernel of a. That is, the complex ^ Gi ^° '""^^'^^ q ^° ''"°'^) p jg exact. 

In both examples, there exist isolated excitations. In the toric code, the isolated excitation 
can be (topologically) nontrivial since the electric charge is not in ime. On the contrary, in 
2D Ising model, any isolated excitation is actually created by an operator of finite support 
because any excitation created by some Pauh operator appears as several connected loops. 
This difference motivates the following definition for charges. 

Let P be the module of Pauli operators of possibly infinite support. Similarly, let E be the 
module of virtual excitations of possibly infinitely many terms. Formally, P is the module 
of all 2g-tuples of functions on the translation- group, and E is that of all t-tuples. Clearly, 
pep and E O E. The containment is strict if and only if the translation-group is infinite. 
Since e consists of Laurent polynomials with finitely many terms, e : P ^ E extends to a 
map from P to E. 

Definition 9. A (topological) charge e — e{p) e is an excitation of finite energy (an 
element of the virtual excitation module) created by a Pauli operator p e P of possibly 
infinite support. A charge e is called trivial if e e e(P). 

By definition, the set of all charges modulo trivial ones is in one-to-one correspondence with 
the superselection sectors. According to the definition, any charge of 2D Ising model is 
trivial. A nontrivial charge may appear due to the following fractal generators. 

Definition 10. We call zero-divisors on cokere as fractal generators. In other words, for a 
complex G P E, eai element / e P\{0} is a fractal generator ii there exists v e £^\ime 
such that fv G ime. 

There is a natural reason the fractal generator deserves its name. Consider a code Hamil- 
tonian with a single type of interaction: t = 1. So each configuration of excitations is 
described by one Laurent polynomial. For example, in two dimensions, f — 1 + x + y = e{p) 
represents three excitations, one at the origin of the lattice and the others at (1,0) and 
(0, 1) created by a Pauli operator represented by p. In order to avoid repeating phrase, let 
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US call each element of the Pauli module a Pauli operator, and instead of using multiplica- 
tive notation we use module operation + to mean the product of the corresponding Pauli 
operators. 

Consider the Pauli operator fp = p + xp + yp&P. It describes the Pauli operator p at the 
origin multiplied by the translations of p at (1, 0) and at (0, 1). So fp consists of three copies 
of p. This Pauli operator maps the ground state to the excited state p = 1 + + y^. The 
number of excitations is still three, but the excitations at (1, 0), (0, 1) have been replaced by 
those at (2, 0), (0, 2). Similarly, the Pauli operator P^^p = P{fp) consists of three copies of 
fp, or 3^ copies of p. The excited state created by pp is = {PY = 1 + + y"^ ■ Still it 
has three excitations, but they are further apart. The Pauli operator consists of 3" 

copies of p in a self-similar way, and the excited state caused by P"~^p consists of a constant 
number of excitations. More generally, if there are t > 1 types of terms in the Hamiltonian, 
the excitations are described by a t x 1 matrix. If it happens to be of form fv for some 
1 7^ / G -R, there is a family of Pauli operators f'^"~^p with self-similar support such that 
it only creates a bounded number of excitations. An obvious but uninteresting way to have 
such a situation is to put fv = e{fp') for a Pauli operator p' where v = e{p'). Our definition 
avoids this triviality by requiring f ^ ime. 

Proposition 5.1. Suppose cokere 7^ 0. Then, the following are equivalent: 

• There does not exist a fractal generator. 

• cokere is torsion-free. 

• There exists a free R-module E' of finite rank such that 

P ^ E' 

is exact. 

Proof. The first two are equivalent by definition. The sequence above is exact if and only if 
— 7- cokere — )■ E' is exact. Since cokere has a finite free resolution, the second is equivalent 
to the third [231 16.33]. □ 

The following theorem states that the fractal operators produces all nontrivial charges. 

Theorem 1. Suppose A = is the translation- group of the underlying lattice. The set of 
all charges modulo trivial ones is in one-to-one correspondence with the torsion submodule 
of coker e. 

Proof. For a module M, let T{M) denote the torsion submodule of M: 

T{M) = {m G M I 3 r e \ {0} such that rm = 0} 

Suppose first that T(cokere) = 0. We claim that in this case there is no nontrivial charge. 
Let e = e{p) G -E be a charge, where p & P. By Proposition 15. II we have an exact sequence 
of finitely generated free modules P E ^ Ei. Since the entries of the matrix ei are finite 
linear combinations of translation-group elements, the complex extends to the complex of 
modules of tuples of functions on the translation-group. 

P ^ E ^ El 

(This extension may not be exact.) Then, ei(e) = ei(e(p)) = since ei o e = identically. 
But, e G -E, and therefore, e G kerei r\ E = t{P). It means that e is a trivial charge, i.e., e 
maps to zero in coker e, and proves the claim. 
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Now, allow coker(P E) to contain torsion elements. Q = (coker e) /T(coker e) is torsion- 
free, and is finitely presented as Q = coker (e' : P' E) where P' is a finitely generated 
free module. In fact, we may choose e' by adding more columns representing the generators 
of the torsion submodule of coker e to the matrix e. Then, P can be regarded as a direct 
summand of P'. 

Let e = e{p) & E he any charge. Since the matrix e' contains e as submatrix, we may write 
e = e'{p) G E. Since T(cokere') = 0, we see by the first part of the proof that e = e'(p') 
for some p' G P'. Then, e maps to zero in Q, and it follows that e maps into T(cokere) in 
coker e. In other words, the equivalence class of e modulo trivial charges is a torsion element 
of coker e. 

Conversely, we have to prove that for every element e & E such that fe = e{p) for some 
f E R\ {0} and p E P, there exists p E P such that e = e(p). Here, P is the module of all 
2g-tuples of F2-valued functions on the translation-group. Consider the lexicographic total 
order on Z^. It induces a total order on the monomials of R. Choose the smallest term /o 
of /. By multiplying /q"^, we may assume /o = 1. (If D = 1, f would be a polynomial of 
nonnegative exponents with the lowest order term 1.) 

We claim that the sequence 



(4) 



converges to / G -R, where R is the set of all F2-valued functions on A. Given the claim, 
since p"e = e + (/ — l)^"e = e(/^" ■ ■ ■ Pfp) where p E P, we conclude that e = e{fp) is a 
charge. 

If / is of nonnegative exponents, and hence / G S = ¥2[xi, . . . ,xd], then the claim is 
clearly true. Indeed, the positive degree terms of = ! + (/ — 1)^" are in the ideal power 
{xi, . . . yXDy C S. Therefore, the sequence Eq. (jl]) converges in the formal power series 
ring F2[[xi, . . . ,xd]], which can be regarded as a subset of R. If / is not of nonnegative 
exponents, one can introduce the following change of basis of the lattice such that / 
becomes of nonnegative exponents. In other words, the sequence Eq. (j3]) is in fact contained 
in a ring that is isomorphic to the formal power series ring where the convergence is clear. 

For any nonnegative integers mi, . . . ,mo-i, define a linear transformation 



Cm C(' 



mi,m2,...,mo-i) 



\aD/ 



/a[\ 



\a'J 



1 1 








••• 0\ 
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mi 


m2 


1 





\mi 


7712 




mo-i 1/ 



(ai\ 
\aD/ 



on Z 



D 



Cm induces the map x"^ ■ • ■ i— x^'^ ■ ■ ■ Xj^ on R. Let u = x"^ ■ ■ ■ be an arbitrary term 
of / other than 1, so n > L That is, for the smallest i G {1, . . . ,D} such that aj 7^ 0, one 
has ttj > due to the lexicographic order. Hence, if we choose rrij large enough and set 
= (j 7^ then Cm(^) has nonnegative exponents. Since any Cm maps a nonnegative 
exponent term to a nonnegative exponent term, and there are only finitely many terms in 
/, it follows that there is a finite composition ( of Cm's which maps / to a polynomial of 
nonnegative exponents. □ 

Since a nontrivial charge v has finite size anyway (the maximum exponent minus the 
minimum exponent of the Laurent polynomials in the t x 1 matrix v), we can say that the 
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charge v is point-like. Moreover, we shall have a description how the point-like charge can 
be separated from the other by a local process. By the local process we mean a sequence of 
Pauli operators [[oi, . . . , o„]] such that Oj+i — Oj is a monomial. The number of excitations, 
i.e., energy, at an instant i will be the number of terms in e(oj). 

Theorem 2. // there is a fractal generator of a code Hamiltonian, then for all sufficiently 
large r, there is a local process starting from the identity by which a point-like charge is 
separated from the other excitations by distance at least T . One can choose the local process 
in such a way that at any intermediate step there are at most or excitations for some constant 
c independent of r. 

For notational simplicity, we denote the local process [[oi, . . . , o„]] by 

S = [Oi, 02 - Oi, 03 - 02,..., On- 0„_i]. 

It is a recipe to construct o„, consisted of single qubit operators. o„ can be expressed as 
"o„ = / s" , the sum of all elements in the recipe. 

Proof. Let / be a fractal generator, and put fv = e{p) where f ^ ime. We already know v 
is a point-like nontrivial charge. Write 

n I 
1=1 i=l 

where each of Pi and fi is a monomial. Let Sq = [0,pi,p2, . . . be a recipe for constructing 
P', J So = p. Given Sj, define inductively 

Si+l = ifl ■ Si) O (/f ■ Si) 0---0 (/f ■ Si) 

where o denotes the concatenation and fi ■ [ui, . . . = [fiUi, . . . , It is clear that 

Si+i constructs the Pauli operator 

J sr = f''-' J 3.-1 = f'^-'f'^" J Sr-2 = f^-'^-^-'^-^^ J So = f'^'p 

whose image under e is f'^'^v. Thus, if r is large enough so that 2^ is greater than the size of 
V, the configuration of excitations is precisely / copies of v. The distance between f 's is at 
least 2^ minus twice the size of v. 

Therefore, there is a constant e > such that for any r > the energy of f'^'^v G -E is < e. 
Let A(r) be the maximum energy during the process s^. We prove by induction on r that 

A(r) < el{r + 1). 

When r = 0, it is trivial. In s.+i, the energy is < A(r) until ff'sr is finished. At the end 
of ff Sr, the energy is < e. During the subsequent /I'^s., the energy is < A(r) + e, and at 
the end of {fY Sr) o (/I'^Sr), the energy is < 2e. During the subsequent ff Sr, the energy is 
< A(r) + je. Therefore, 

A(r + 1) < A(r) + el< el{r + 2) 
by the induction hypothesis. □ 

Note that the notion of fractal generators includes that of 'string operators'. In fact, a 
fractal generator that contains exactly two terms gives a family of nontrivial string segments 
of unbounded length, as defined in [7]. The proof here mimics the explicit construction of 
fractal operators in [8]. 
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Below, we point out a couple of sufficient conditions for nontrivial charges, or equivalently, 
fractal generators to exist. 

Proposition 5.2. For code Hamiltonians, the existence of a fractal generator is a property 
of an equivalence class of Hamiltonians in the sense of Definitionl^ 

Proof. Suppose ima = ima'. Each column of a' is a i?-linear combination of those of a, 
and vice versa. Thus, there is a matrix B and B' such that e' = Be and e = B'e'. BB' 
and B'B are identity on ime' and ime respectively. In particular, B' and B are injective on 
ime' and ime respectively. Suppose / is a fractal generator for e, i.e., fv = ep ^ 0. Then, 

7^ Bfv = fBv = Be{p) = t'{p). If Bv G ime', then v = B'Bv G ime, a contradiction. 
Therefore, / is also a fractal generator for e'. By symmetry, a fractal generator for e' is a 
fractal generator for e, too. 

Suppose i?' C i? is a coarse-grained base ring. If cokere is torsion-free as an i?-module, 
then so it is as an i?'-module. If / G -R is a fractal generator, the determinant of / as a 
matrix over R' is a fractal generator. 

A symplectic transformation or tensor ing ancillas does not change coker e. □ 

Proposition 5.3. For any ring S and t > I, ifO^S^^ S^* 5* is exact and ^ S, 
coker is not torsion-free. In particular, for a degenerate exact code Hamiltonian, if o is 
injective, then there exists a fractal generator. 

Proof. By Proposition 13.21 rank0 = t. Since C It{(f)) C S" is the initial Fitting ideal, we 
have 7^ ann coker 7^ S. That is, coker is not torsion-free. 

For the second statement, set S = R. If a is injective, we have an exact sequence 

O^G^P^E. 

By Remark [21 t = rank G = rank a = rank e = q. □ 

Proposition 5.4. Suppose the characteristic dimension is D — 2 for a degenerate exact code 
Hamiltonian. Then, there exists a fractal generator. 

Proof. Suppose on the contrary there are no fractal generators. Then, by Proposition 15. 

G^P^E^E' 

is exact for some finitely generated free module E'. Since coker e has finite free resolution by 
Lemma \7.1\ Proposition 13.21 implies codim I{a) > 3 unless I{a) = R. But, codim/(cr) = 2 
and I{a) ^ Rhj Corollary 14.21 This is a contradiction. □ 

6. One dimension 

The group algebra R = F2[x,x] for the one dimensional lattice Z is a Euclidean domain 
where the degree of a polynomial is defined to be the maximum exponent minus the minimum 
exponent. (In particular, any monomial has degree 0.) Given two polynomials /, g in R, one 
can find their gcd by the Euclid's algorithm. It can be viewed as a column operation on the 
1x2 matrix (/ g^ . Similarly, one can find gcd of n polynomials by column operations on 

1 X n matrix 

(/i /2 ■ ■ ■ fn) ■ 
The resulting matrix after the Euclid's algorithm will be 

(gcd(/i,...,/„) ■■■ 0). 
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Given a matrix M of univariate polynomials, we can apply Euclid's algorithm to the first 
row and first column by elementary row and column operations in such a way that the degree 
of (1, l)-cntry Mn decreases unless all other entries in the first row and column are divisible 
by Mil. Since the degree cannot decrease forever, this process must end with all entries in 
the first row and column being zero except Mn. By induction on the number of rows or 
columns, we conclude that M can be transformed to a diagonal matrix by the elementary 
row and column operations. This is known as the Smith's algorithm. 

The following is a consequence of the finiteness of the ground field. 

/(^) >^ 

Lemma 6.1. Let ¥ be a finite field and S = ¥[x] be a polynomial ring. Letcj): S S 
be a 1 X 1 matrix such that /(O) 0. (f) can be viewed as an n x n matrix acting on the 
free S'-module S where S' — ¥[x'] and x' — x^ . Then, for some n > 1, the matrix (f) is 
transformed by elementary row and column operations into a diagonal matrix with entries 1 
or x' — 1. The number of x' — 1 entries in the transformed 4> is equal to the degree of f. 

Proof. The sphtting field F of f{x) is a finite extension of F. Since F is finite, every root of 
f{x) is a root of — 1 for some n' > 1. Choose an integer p > 1 such that 2^ is greater 
than any multiplicity of the roots of f{x). Then, clearly f{x) divides {x"" — 1)"^^ = x'^''"' — 1. 
Let n be the smallest positive integer such that f{x) divides a;" — 1. 

Consider the coarse- graining by S' — ¥[x'] where x' — x". 5" is a free 5"- module of rank 
n. and (/) is now an endomorphism of the module S represented as an n x n matrix. Since 
f{x)g{x) — x^ — 1 for some g[x) e F[x], we have 

AB = {x' - l)idn 

where x' = x", and A, B are the matrix representation of f{x) and g{x) respectively as 
endomorphisms. A and B have polynomial entries in variable x' . The determinants of ^4, 
are nonzero for their product is {x' — 1)" ^ 0. Let Ei and E2 be the products of elementary 
matrices such that A' — E1AE2 is diagonal. Such matrices exist by the Smith's algorithm. 
Put B' = E^^BEi^. Then, 

A'B' = EiAE2E2^BE{^ = E^ABE^^ = {x' - l)id„. 

Since A' and In are diagonal of non-vanishing entries, B' must be diagonal, too. It follows 
that the diagonal entries of A' divides {x' — 1); that is, they are 1 or x' — 1. 

The number of entries x' — 1 can be counted by considering S/ {f{x)) as an F- vector space. 
It is clear that diuif S / {f (x)) = degf{x). S/{f{x)) — coker0 viewed as a 5"-module is 
isomorphic to 5""/imA', the vector space dimension of which is precisely the number of 
x' — 1 entries in □ 

Theorem 3. If K='L, any system governed by a code Hamiltonian is equivalent to finitely 
many copies of Ising models, plus some non-interacting qubits. In particular, the topological 
order condition is never satisfied. 

We will make use of the elementary symplectic transformations and coarse-graining to 
deform o" to a familiar form. Recall that for any elementary row-addition E on the upper 
block of (7 there is a unique symplectic transformation that restricts to E. 
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Proof. Applying Smith's algorithm to the first row and the first column oi 2q x t matrix a, 
one gets 

//i o\ 
A 



91 



\ 



92 



by elementary symplectic transformations. Let 1 < z < j < g be integers. If some + 
entry is not divisible by /i, apply Hadamard on j-th qubit to bring (g + j)-th row to the 
upper block, and then run Euclid's algorithm again to reduce the degree of (1, l)-entry. The 
degree is a positive integer, so this process must end after a finite number of iteration. Now 
every {q + j, l)-entry is divisible by /i and hence can be made to be by the controUed- 
NOT-Hadamard: 

//i o\ 

A 



9i 



92 



Further we may assume deg /i < deg gi. Since a'^Xqa = 0, we have a commutativity condition 



fi9i - 9ifi = 0. 

Write /i = ax"" + ■ ■ ■ + (3x^ and gi = •yx'^ + ■ ■ ■ + Sx'^ where a < b and c < d and a, f3,'y,S ^ 0. 
Then, figi = f3'yx'^~^ + ■ ■ ■ + a6x'^~°'. Since figi = figi, it must hold that — (c — 6) = d — a 
and a6 = f3'y. Since deg/i < deg^^i, we have d — b = — (c — a) > 0. The controlled-Phase 
Ei+q^i{—{x'''~^ + x'^~°')5 / (3) will decrease the degree of gi by two, which eventually becomes 
smaller than deg/i. One may then apply Hadamard to swap /i and gi. Since the degree of 
(1, l)-entry cannot decrease forever, the process must end with gi = 0. 

The commutativity condition between i-th{i > 1) column and the first is ficji = 0. Since 
fi ^ 0, we get gi = 0: 

//i o\ 

A 







Continuing, we transform a into a diagonal matrix. (We have shown that a can be trans- 
formed via elementary symplectic transformations to the Smith normal form.) 

Now the Hamiltonian is a sum of non-interacting purely classical spin chains plus some 
non-interacting qubits {fi = 0). It remains to classify classical spin chains whose stabilizer 
module is generated by 

(/) 

where we omitted the lower half block. We can always choose / = f{x) such that f{x) has 
only non-negative exponents and /(O) 7^ since x is a unit in R. Lemma [6.11 says that (/) 
becomes a diagonal matrix of entries 1 or x' — 1 after a suitable coarse-graining followed 
by a symplectic transformation and column operations. 1 describes the ancilla qubits, and 
x' — 1 = x' + 1 does the Ising model. □ 

According to the proof, for the interaction range w of classical Hamiltonian, the amount 
of coarse-graining n needed to map it to the Ising models can be exponential in w. One 
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cannot naively reduce n; if n = 2^" — 1, the minimal polynomial of a primitive n-th root of 
unity has degree w. 

7. Two DIMENSIONS 

In the following two sections we will be mainly interested in exact code Hamiltonians. If 
D = 2, the lattice is A = Z^, and our base ring is R = ¥2[x, x, y, y]. 

The following asserts that the local relations — a few terms in the Hamiltonian that 
multiply to identity in a nontrivial way as in 2D Ising model, or the kernel of a — among 
the terms in a code Hamiltonian, can be completely removed for exact Hamiltonians in two 
dimensions. We prove a more general version. 

Lemma 7.1. If G P E is exact over R = F2[xi, Xi, . . . , x^, x^], There exists a' : 
G' ^ P such that im a' = ima and 

^ Gd-2 y Gi ^ G' ^ P ^ E 

is an exact sequence of free R-modules. If D = 2, one can choose o' to he infective. 

We make use of a constructive version of Hilbert's syzygy theorem via Grobner basis. 

Proposition 7.2. [I9l Theorem 15.10, Corollary 15.11] Let {^fi, . . . he a Grohner hasis 
of a suhmodule of a free module Mq over a polynomial ring. Then, the S-polynomials Tij 
of {di} in the free module Mi = ^^^^Sci generate the syzygies for {gi}. If the variahle 
Xi, . . . ,Xs are ahsent from the initial terms of gi, one can define a monomial order on Mi 
such that hsent from the initial terms of Tij. If all variahles are ahsent from 

the initial terms of gi, then Mq/ {gi, . . . , (?„) is free. 



Proof of Lemma 7. 1 . Without loss of generality assume that the tx2q matrix e have entries 
with nonnegative exponents so they are in fact polynomials, not Laurent polynomials. Below, 
every module is over the polynomial ring S = F2[xi, . . . ,xd] unless otherwise noted. Let E+ 
be the free 5- module of rank equal to rank/j E. 

If gi,--- ,g2q are the columns of e, apply Buchberger's algorithm to obtain a Grobner 
basis gi, - ■ ■ , g2q, ■ ■ ■ ,gn of ime. Let e' be the matrix whose columns are gi, . . . , gn- We 
regard e' as a map Mq — )■ By Proposition 17. 2[ the initial terms of the syzygy generators 
(S-polynomials) Tij for {gi} lacks the variable xi. Writing each Tij in a column of a matrix 
Tl, we have a map ri : Mi — t- Mq. 

By induction on we have an exact sequence 

Md ^ Md-1 ■ ■ ■ ^ Mo 4 

of free 5'-modules, where the initial terms of columns of t^ lack all the variables. By Propo- 
sition [7[2] again, M']j_i = Md-i/ iuiTT-^ is free. Since kerTo-i = imr/j, we have 

^ M'^^i • • ■ ^ Mo 4 

Since g2q+i, ■ ■ ■ ,gn are S'-linear combinations of gi, . . . , g2q, there is a basis change of Mo so 
that the matrix representation of e' becomes 

e' ^ (e 0) . 

With respect to this basis of Mo, the matrix of ti is 



n = 



rid 
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where tiu is the upper 2q x t' submatrix. Since kere' = imri, The first row r of Tid should 
generate 1 G 5. (This property is called unimodularity.) Quillen-Suslin theorem [211 Chap- 
ter XXI Theorem 3.5] states that there exists a basis change of Mi such that r becomes 
(l ■ ■ ■ O). Then, by some basis change of Mq, one can make 



where is a submatrix. By induction on the number of rows in Tid, we deduce that the 
matrix of ti can be brought to 



e' = (e 0) , n 







Note that ecr" = and ea' = 0. The basis change of Mq hj [ ^ ^ ) gives 



-a 



The kernel of J determines kerri = imr2. Let M[ denote the projection of Mi such that 
the sequence 

^ M'j^_^ >M2^ M[ ^ M'^ 4 E+ 

of free S'-modules is exact. 

Taking the ring of fractions with respect to the multiplicatively closed set 

t/ = K---a;*^|zi,...,z^>0}, 

we finally obtain the desired exact sequence over U~^S = R with P = U^^M'q and E = 
U~^E^. Since imo" = kere, we have imcr' = imo". □ 

Lemma 7.3. Let R be a Laurent polynomial ring in D variables over a finite field ¥, and N 
be a module over R. Suppose J = annfiN is a proper ideal such that dimR/J = 0. Then, 
there exists an integer L > 1 such that 

anuN = (xf - 1, . . . , - 1) C i?' 

where R' = F[xf', , . . . , xj^, x^] is a subring of R. 

This is a variant of Lemma 16.11 

Proof. Since i? is a finitely generated algebra over a field, for any maximal ideal m of i?, 
the field R/xn is a finite extension of F (NuUstellensatz [191 Theorem 4.19]). Hence, R/m 
is a finite field. Since Xi is a unit in R, the image G -R/m of Xi is nonzero. Oj being an 
element of finite field, a power of is 1. Therefore, there is a positive integer n such that 

= {^i — 1, . . . , — 1) Cm. Since x"" — 1 divides x""" — 1, we see that there exists n>l 
such that b„ C rrii n m2 for any two maximal ideals rrii,m2. One extends this by induction 
to any finite number of maximal ideals. 

Since dim R/ J = 0, any prime ideal of R/ J is maximal and the Artinian ring R/ J has only 
finitely many maximal ideals, rad J is then the intersection of the contractions (pull-backs) 
of these finitely many maximal ideals. Therefore, there is n > 1 such that 

br, C rad J. 
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Since R is Noetherian, (rad J)^"^ C J for some r > where p is the characteristic of F. 
Hence, we have 

Kpr C C (rad J)P' C J. 

Let L = nf. If R' = ¥[ Xi, Xi, . . . , xj^, ann/j/ A/" is nothing but J Pi i?'. We have just 
shown bL^R' C JnR'. Since J is a proper ideal, we have 1 ^ JnR'. Thus, bifli?' = Jfli?' 
since b/, fl i?' is maximal in R'. □ 

Theorem 4. For any two dimensional degenerate exact code Hamiltonian, there exists an 
equivalent Hamiltonian such that 

anncoker e = (x — 1, ?/ — 1). 

Thus, cokere is a torsion module. 

Proof. By Lemma 17.11 we can find an equivalent Hamiltonian such that the generating map 
a for its stabilizer module is injective: 

^ G ^ P. 

Let t be the rank of G. The exactness condition says 

is exact where e = 0"^^ and E has rank t. Applying Proposition 13. 2^ since /(a) = /(e) 
and hence in particular codim/(o") = codim/(e), we have that q = t and codim/(e) > 2 if 
J(e) ^ R. But, /(e) ^ P by Corollary i^l 

Since q = t, J(e) is equal to the initial Fitting ideal, and therefore has the same radical 
as the annihilator of cokere = E/ime. (See |19i Proposition 20.7] or [211 Chapter XIX 
Proposition 2.5].) In particular, dim i?/(ann cokere) = 0. Apply Lemma [7^ to conclude the 
proof. □ 

An interpretation of the theorem is the following. 

For systems of qubits. Theorem H] says that x + 1 and y + I are in ann coker e. In other 
words, any element f of is a charge, and a pair of f 's of distance 1 apart can be created by 
a local operator. Equivalently, v can be translated by distance 1 by the local operator. Since 
translation by distance 1 generates all translations of the lattice, we see that any excitation 
can be moved through the system by some sequence of local operators. This is exactly what 
happens in the 2D toric code: Any excited state is described by a configuration of magnetic 
and electric charge, which can be moved to a different position by a string operator. 

Moreover, since {x — l,y — 1) = ann cokere, the action oi x,y G i? on cokere is the 
same as the identity action. Therefore, the i?-module coker e is completely determined up to 
isomorphism by its dimension k as an F2-vector space. The module K{L) of Pauli operators 
acting on the ground space (logical operators), can be viewed as K{L) = Tori (coker e, R/bi). 
Thus, K{L) is determined by k up to isomorphisms. This implies that the translations of 
a logical operator are all equivalent. It is not too obvious at this moment whether the 
symplectic structure, or the commutation relations among the logical operators, of K{L) is 
also completely determined. 

A similar result is reported in flUi [TT] . 
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8. Three dimensions 

In the previous section, we derived a consequence of the exactness of code Hamiltonians. 
The two-dimensional Hamiltonian was special so we were able to characterize the behavior 
of the charges more or less completely. Here, we prove a weaker property of three dimensions 
that there must exist a nontrivial charge for any exact code Hamiltonian. It follows from 
Theorems [1] and [2] that such a charge can spread through the system by surmounting the 
logarithmic energy barrier. 



Lemma 8.1. Suppose D = 3, 

O^Gi^G^P E 
is exact, and I (a) C m = (x — 1, y — 1, 2 — 1). Then, coker e is not torsion-free. 
Proof. Suppose on the contrary coker e is torsion-free. We have an exact sequence 

If Gi =0, Proposition 15.31 implies the conclusion. So we assume Gi ^ 0, and therefore we 
have /(cTi) = R. 

Let us localize the sequence at m to have /((Ti)^ = Rm- Since rank(G'i)m = rank(cri)n^, the 
matrix of (ai)m becomes 



for some basis of {Gi)m and Gm- See the proof of Lemma 14. 1[ In other words, there is an 
invertible matrix B G GLfxt(-Rm) such that 

a^B = (a 0) 

where a is the 2q x t' submatrix. Note that the antipode map is a well-defined automorphism 
of Rra since m = m. 

Since e = cr^Xq, we have 



(5) = A 

Therefore, we get a new exact sequence 



m 



where G' = G^/ im((Ji)n, is a free i^m-niodule and t' = rankC. It is clear that ranke = 
rank 5". Setting S = Rm in Proposition 15.31 implies that coker e is not torsion-free. But, since 
we are assuming coker is torsion- free, coker a''" is also torsion- free by Eq. (|5]). This is a 
contradiction. □ 

Theorem 5. For any three-dimensional, degenerate and locally topologically ordered code 
Hamiltonian, there exists a fractal generator. 

Proof. By Lemma 17. H there exists an equivalent Hamiltonian such that 

^ Gi ^ G A P E 
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is exact. The existence of a fractal generator is a property of the equivalence class by 
Proposition 15.21 If we can find a coarse-graining such that /(cr') C [x' — l,y' — l,z' — 1), 
then Lemma 18.11 shall imply the conclusion. 

Recall that and o"/, denote the induced maps by factoring out = — y^ — l,z'" — l). 
See Sec. HI There exists L such that K{L) = kere^/ imcxi 7^ by Corollary 14.21 Consider 



the coarse-grain by x' 



x^, y' 



. Let R' = F2[x''^^ , y'"^^ , z'"^^] denote the coarse- 



grained base ring. If K'{L') denotes kere^,/ imo"^, as -R'-module, we see that -^''(1) = K{L) 
as F2-vector space. In particular, K'{1) ^ 0. Put m = (x' — l,y' — 1, — 1) = C R' . 
Then, K'{1)^ = K'{1) ^ 0. By Lemma ill we have J(a') Cm. □ 



9. More examples 

Example 4 (Toric codes in higher dimensions). Any higher dimensional toric code can be 
treated similarly as for two dimensional case. In three dimensions one associates each site 
with q = 3 qubits. It is easily checked that 



3D-toric 



1 + y 
1 + z 





















V 











l + z 
1 + y 



l + z 


1 + x 



1+y 
1 + x 




Both two- and three-dimensional toric codes have the property that coker e is not torsion- 
free. However, in two dimensions any element of is a physical charge, whereas in three 
dimensions E contains physically irrelevant elements. Note that in both cases, 1 + x and 1 + y 
are fractal generators. Being consisted of two terms, they generate the 'string operators'. 

The 4D toric code [1] has a^-type interaction and a^-type interaction. Originally the 
qubits are placed on every plaquette of 4D hypercubic lattice; instead we place q = 6 qubits 
on each site. The generating map a for the stabilizer module is written as a 12 x 8-matrix 
{t = 8) 



C"4D-toric 



0-x 









where 





1 + x 








\ 


1 + w 








1 + 


X 


l + z 





1 + x 










l + z 


1 + y 










1 + w 





1+y 


V 





1 + w 


1 + 




( ° 





1 + w 


1 + 







l + z 


1+y 










1 + w 





1 + 


y 


1 + w 








1 + 


X 


l + z 





1 + x 







Vi + y 


1 + x 








) 



Note the bar on oz- 
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Theorem [5] does not prevent the absence of a fractal generator in four or higher dimensions. 
Indeed, this 4D toric code lacks any fractal generator. To see this, it is enough to consider 

oz since coker a 



X 



coker 0"^ as -R-modules. If 



= (1 + X l + y l + z 1 + w) : R"^ ^ R, 



then 



R 



R 



is exact. (A direct way to check it is to compute S-polynomials [IHl Chapter 15] of the entries 
of ei, and to verify that they all are in the rows of cr^.) Hence, coker is torsion-free by 
Proposition 15.11 

For the toric codes in any dimensions, a has nonzero entries of form Xi — \. The radical of 
the associated ideal J(cr) is equal to m = (a;i — 1, . . . , — 1). So m is the only maximal ideal 
of R that contains /(cr). The characteristic dimension is zero. If 2 f L, since {^l)^^ = TTim, 
{<^L)m is a zero matrix. Any other localization of ai does not contribute to dimpj K{L) by 
Lemma 14. 1[ Therefore, if 2 f L, K{L) has constant vector space dimension independent of 
L; the ground state degeneracy is independent of system size. 

Example 5 (Wen plaquette [25]). This model consists of a single type of interaction (t = 
q=l) 



X-Y 
I I 
Y-X 



C"We: 



1 + X + y + xy 
1 + xy 



where X,Y are abbreviations of (Tx,(Ty. It is known to be equivalent to the 2D toric code. 
Take the coarse-graining given by R' = ¥2[x', y', x', y'] where 



X 



xy, 



y' = y'. 



(The coarse-graining considered in this example is intended to demonstrate a non-square 
blocking of the old lattice to obtain a 'tilted' new lattice, and is by no means special.) As 
an i?'-module, R is free with basis {1, y}. With the identification R = {R' ■ 1) © {R' ■ y), we 
have X ■ 1 = x' ■ y, X ■ y = x'y' ■ 1, and y ■ 1 = 1 ■ y, y ■ y = y' ■ 1. Hence, x and y act on 
i?'-modules as the matrix- multiplications on the left: 



X 



x'y' 



y 



y' 

1 



Identifying 



R- = [(R' ■ 1) © (i?' ■ y)] © ■ ■ ■ © [{R' ■ 1] 

our new a on the coarse-grained lattice becomes 

/l + x'y' y' + x'y'\ 
1 + x' 1 + x'y 



{R' 



a 



1 + x'y' 
\ 1+x'y'J 
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/ ^2,4(1)^ 
El,3(l) 



/ y' + x'y'\ 
1 + x' 
1 + x'y' 

y 1 + x'y'J 



E3.2{y') 



+ x'y'\ 



( 



col. 2 



1 +x' 



1 + x' 
l + y' 

1 + ^7 



1^3 



v 



( 

1 + 

1 + 

y : x'y' + x'y 

\ 



/I 

1 + x' 

1 + x' 

V 1+ry 



which is exactly the 2D toric code. 







Example 6 (Chamon model [261 ETj). This three-dimensional model consists of single type 
of term in the Hamiltonian. The generating map is 



C"Cha 



X + X + ?/ + ?/ 

z^z+y^y 



Since 



atAi = (1 + xy) 





X + ?/ 



1 + xy is a fractal generator. Consisted of two terms, it generates a string operator. The 
degeneracy can be calculated using Corollary 14.51 Assume all the three linear dimensions of 
the system are even. Put 



S = R/{x + x + y + y,z + z + y + y,x''-l,y'"'-l,z'''-l). 

Then, the loga of the degeneracy is k = dim^^ S. In S, we have x + x = y + y = z + z. Since 
S has characteristic 2, it holds that 



{w + w-^){wP + wP-^ + ■■■ + w-P) 



for p > 1 and w = x,y, z. By induction on p, we see that + w ^ is a polynomial in 
w + w~^. Therefore, 

xP + = yP + f = z^ + zP 

for all p > 1 in S. Put g = gcd(/, m, n). Since x' + x^' = ?/'" + y~"^ = 2;" + = in S", we 
have x^ + x~3 = + = + z'^ = 0. 

Applying Buchberger's criterion with respect to the lexicographic order in which x < y < z, 
we see that 

S = F2[x, y, z\l{z^ + zx^^~^ + 2x + 1, + yx^^~^ + yx + 1, x^^ + 1) 
is expressed with a Grobner basis. Therefore, 

k = dim]F2 S = 8 gcd(Z, m, n). 
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Example 7 (Cubic Code). The Hamiltonian of code 1 in [7] is the translation-invariant 
negative sum of the following two types of interaction terms: 



IZ 



Zl 



IZ 



ZZ 



u 



Zl 



IZ 



IX 



XI 



II 



XX 



XI 



IX 



xz 



Zl 



IX 



XI 



xyz 



yz 



X 



xy 



^cubic—code 



\ 



Here, the third cube specifies the coordinate system of the simple cubic lattice. The corre- 
sponding generating map for the stabilizer module is 

fl + xy + yz + zx \ 

1 + X + y + z 

1 + X + y + z 

1 + xy + yz + zx ^ 

The associated ideal is contained in a prime ideal of codimension 2: 

-^(cr) ^ (1 + a; + ?/ + 2;, 1 + a;?/ + + zx) = p. 

Since codim /(a) > 2, the characteristic dimension is 1. Since coker ecubic-code = -R/p © -R/p, 
any nonzero element of p is a fractal generator. 

Let us explicitly calculate the ground state degeneracy when the Hamiltonian is defined 
on L X L X L cubic lattice with periodic boundary conditions. By Corollary 14. 5[ 

k = dimF2 R/ip + bi) © R/{p + bi) = 2 dim^^ R/{p + bi). 

So the calculation of ground state degeneracy comes down to the calculation of 

d = dimp^ T'/p 

where T' = F2[x, y, zj/ix""^ - 1, y"^ _ 1^ z""^ _ 1). 

We may extend the scalar field to any extension field without changing d. Let F be the 
algebraic closure of F2 and let 

T = F[a;, y, - 1, y"^ - 1, z"^ - 1) 

be an Artinian ring. By Proposition 14.31 it suffices to calculate for each maximal ideal m of 
T the vector space dimension 

dm = dimF(T/p)n, 
of the localized rings, and sum them up. 

Suppose rii, 71,2,^3 > 1- By Nullstellensatz, any maximal ideal of T is of form m = 
(x — xo,y — yo, z — zo) where Xq^ = y^^ = Zq^ = 1. (If rii = n2 = = 1, then T becomes a 
field, and there is no maximal ideal other than zero.) Put rii = 2''n- where n'^ is not divisible 
by 2. Since the polynomial x^^ — 1 contains the factor x — xq with multiplicity 2'^, it follows 
that 



T 



¥[x,y,z]J{x''' +a', y'"'+b', z''' + c') 



where a' 



yl \c' = zl \ Hence, (T/p)^ = F[x, z]/I' where 



/' = {x + y + z + l,xy + xz + yz + 1, 
If r = ¥[x,y,z], then d^ = 0. 



X 



2'1 



+ a , y +b , z + c) 
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Without loss of generality, we assume that h < h ^ h- By powering the first two 
generators of we see that {xq, i/o, Zq) must be a solution of them in order for I' not to be 
a unit ideal. Eliminating z and shifting x ^ x + 1, y— >-y + l, our objective is to calculate 
the Grobner basis for the proper ideal 

I ^{x'^ + xy + y^, x^'' + a, y^'' + b) 

where a — a' + 1 and b — b' + 1. So 

dm = dimiFF[a;,y]//. 

One can easily deduce by induction that t/^"" + x'^"'~^{'mx + y) & I for any integer m > 0. 
And b — o'o^'^ for a primitive third root of unity u. So we arrive at 

/ = (y2 _^ + a;^ i/x^''"^ + 6(1 + kco^), x^'^ + a) 

We apply the Buchberger criterion. If a 7^ 0, i.e., xq 7^ 1, then 6 7^ and I = {x + (cu^ + 
l2)y,x^'^ +a), so dm = 2'i 

If a = 6 = 0, then I — (y'^ + yx + x^, yx"^^'^, x^^). The three generators form Grobner 
basis if I2 — h- Thus, in this case, dm = 2^^+^ — 1. If I2 > h, then d^ — 2'^+^. 

To summarize, except for the special point (1, 1, 1) e of the affine space, each point in 
the algebraic set 

X + y + z + 1 = xy + xz + yz + 1 = 
x< - 1 = y< - 1 = - 1 = 

contribute 2'^ to d. The contribution of (1, 1, 1) is either 2^^"'"^ or 2'^"'"^ — 1. The latter occurs 
if and only if /i and the two smallest numbers of factors of 2 in rii, n2, 713, are equal. Let 
do — i^V be the number of points in V. The desired answer is 



V^\{x,y,z)er 



d = 2'i(do-l) + 




1 if/i = /2 
otherwise 



where h < I2 ^ h are the number of factors of 2 in n,. 

The algebraic set defined by {x + y + z + 1, xy + xz + yz + 1) is the union of two isomorphic 
lines intersecting only at x — y — z — 1, one of which is parametrized by x e F as 

(l + x,l+ux,l + oj^x) e F^, 

and another is parametrized as 

(1 + x, 1 + oo'^x, 1 + oox) e F^. 

where a; is a primitive third root of unity. Therefore, the purely geometric number do = 
2di — 1 can be calculated by 

di = deg^ gcd ((1 + x)< + 1,(1 + uxf^ + 1,(1 + uj'^xp + 1) . 

In particular, ii L — rii — n2 — > 1 (including all the factors of 2), 

(6) ^^ = deg^gcd((l + a;)^ + l, (l + a;a;)^ + l, (l + cj^xf + l) 

for any L > 2. The true ground state degeneracy is 2'^ — 4*^. Some special cases are explicitly 
computed: 
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'1 


zfL 


k 




2L 


-1 zfL 


2 




2L 


-5 z/L 










ED- 


Use (a 


+ (3r = 



Corollary 9.1. Let k be log2 of the ground state degeneracy of the cubic code on the cubic 
lattice of size with periodic boundary conditions. Let p > 1 be any integer. Then, 

2P + 1, 
2P, 

4P - 1, 

22P+1 _ i_ 

Proof. Recall k = 2d in Eq. ([6]). Use (a + /3f" = a^" + P"^^ and a;^ ^ w + 1 = 0. □ 

10. Discussion 

There are many natural questions left unanswered. Perhaps, it would be the most in- 
teresting to answer how much the associated ideal I {a) determines about the Hamiltonian. 
Note that the very algebraic set defined by the associated ideal is not invariant under coarse- 
graining. For instance, in the characteristic dimension zero case, the algebraic set can be a 
several points in the affine space, but becomes a single point under a suitable coarse-graining. 

It is reasonable to conceive that the algebraic set is mapped by the affine map (oj) H- (a^) 
under the coarse-graining by x[ = xf. This is true if t = q, so the g-th determinantal ideal 
of e, being the initial Fitting ideal, has the same radical as anncokere. In fact, we have 
implicitly used this idea in the proofs of Lemma 16. H 17. 3[ and Theorem O The case t > q is 
not explicitly handled here. 

Also, it is interesting on its own to prove or disprove that the elementary symplectic 
transformations generate the whole symplectic transformation group. 
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